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ABSTRACT 

We investigate the generation of magnetic field by the Biermann battery in 
cosmological ionization fronts, using new simulations of the reionization of the 
universe by stars in protogalaxies. Two mechanisms are primarily responsible for 
magnetogenesis: i) the breakout of I-fronts from protogalaxies, and ii) the propagation 
of I-fronts through the high density neutral filaments which are part of the cosmic 
web. The first mechanism is dominant prior to overlapping of ionized regions {z ~ 7), 
whereas the second continues to operate even after that epoch. However, after overlap 
the field strength increase is largely due to the gas compression occurring as cosmic 
structures form. As a consequence, the magnetic field at z ~ 5 closely traces the gas 
density, and it is highly ordered on megaparsec scales. The mean mass-weighted field 
strength is Bq ~ 10~^^ G in the simulation box. There is a relatively well-defined, 
nearly linear correlation between Bq and the baryonic mass of virialized objects, with 
Bq « 10-1^ G in the most massive objects (M ^ W^Mq) in our simulations. This is 
a lower limit, as lack of numerical resolution prevents us from following small scale 
dynamical processes which could amplify the field in protogalaxies. Although the field 
strengths we compute are probably adequate as seed fields for a galactic dynamo, 
the field is too small to have had significant effects on galaxy formation, on thermal 
conduction, or on cosmic ray transport in the intergalactic medium. It could, however, 
be observed in the intergalactic medium through innovative methods based on analysis 
of 7-ray burst photon arrival times. 

Subject headings: cosmology: theory - cosmology: large-scale structure of universe - 
galaxies: formation - galaxies: intergalactic medium 



1. Introduction 

The origin of magnetic fields in the universe is one of the most important unsolved problems 
in cosmology, and has possible ramifications for theories of galaxy formation and evolution 
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(Wasserman 1978, Kim et al. 1996, Kulsrud et al. 1997) and for theories of the formation of the 
first stars (Yamada Sz Nishi 1998, Omukai & Nishi 1998, Abel et al. 1998, Nakamura 8z Umemura 
1999, Bromm, Coppi, &: Larson 1999); see also Rees (1987). 

It is widely accepted that astrophysical magnetic fields reached their present state in a 
two stage process; first, the generation of a seed field, and then hydromagnetic processes which 
amplified the field or smoothed it. What are the possible sources for the origin of such seed 
field? The question has been posed repeatedly since the original work by Biermann (1950). 
Generally speaking, at least three different possible scenarios have been put forward. The first one 
assumes that a magnetic field is present ab initio: this does not violate homogeneity but makes 
the cosmological expansion anisotropic, and affects primordial nucleosynthesis via an increase of 
electron density. The field could have been generated either during QCD/electroweak first order 
cosmic phase transitions (Quashnock, Loeb &: Spergel 1989, Vachaspati 1991) or during infiation 
(Turner &; Widrow 1988, Ratra 1992). The upper limits on the present field strength obtained 
from CMB anisotropy measurements (Barrow, Ferreira &; Silk 1997 give -Bq < 7 x 10~^(J7/i^)^/^ G) 
appear to be more stringent than those derived from the Big Bang nucleosynthesis. The latter 
have been recently revised by Grasso & Rubinstein (1995), who found Bq < 3 x 10~^ G. The 
upper limit from CMB measurements comes very close to the minimum value (~ 10^^ G) required 
for magnetic stresses to have been dynamically important (for example in the creation of voids, de 
Araujo & Opher 1997) in the post-recombination epoch. Compression alone, without a dynamo, 
would have brought such fields to the /xG strengths characteristic of magnetic fields in galaxies. It 
is not clear that magnetic fields generated in the early universe possessed the degree of spatial 
coherence seen in galactic magnetic fields (Zweibel & Heiles 1997) but they could have been 
substantially reconfigured by turbulence and differential rotation in galaxies. 

The second type of field generation mechanism would have operated prior to the recombination 
epoch, and was first proposed by Harrison (1970). This process invokes a Biermann battery 
(originally proposed for stars), i.e. nonparallel pressure and density gradients, powered by a 
non-zero vorticity in the primordial fluctuation field. Electrons and ions would have tended to 
spin at different rates due to the drag on electrons by the CMB. An ambipolar electric field would 
have arisen to couple the electrons and ions, and a magnetic field would thus have been created by 
induction. Zeldovich, Ruzmaikin, & Sokoloff (1983) later discussed a variant of this mechanism. 
Even if the vorticity is equivalent to galactic rotation hy z = 10, the magnitude of the resulting 
magnetic field is only of order 10~^^ G. However, vorticity modes are now thought to decay rapidly 
as irrotational modes grow and structure forms; thus this process now appears very unlikely. 

In the third class of schemes, magnetic fields can arise during the epoch of protogalaxy 
formation, through any mechanism which produces a non-potential electric field. It is expected that 
fields produced in this manner must be small, and only act as a seed for a hydromagnetic dynamo. 
Two mechanisms for producing seed fields have been considered. In the first, magnetogenesis 
proceeds primarily in shocks. Davies & Widrow (1999) showed with simple analytic models 
that magnetic fields as high as 10~^^G can be generated during the collapse of protogalaxies. 
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Kulsrud et al. (1997) computed the effect of the Biermann battery in a simulation of large scale 
structure. They found that magnetic fields are generated primarily in shocks associated with the 
gravitational collapse of the structure, with a a characteristic field strength of 10~^^ G. 

In the second mechanism, magnetogenesis takes place in ionization fronts. Subramanian, 
Narasimha, & Chitrc (1994, SNC) pointed out that the Biermann battery would operate in 
cosmological ionization fronts propagating through density irregularities. They made simple 
estimates based on the properties of protogalactic fluctuations, and suggested that fields of a few 
times 10~^° G, with a coherence length of several kpc, could be generated by this mechanism. 
Along similar lines, but exploiting SN explosions rather than ionization fronts, Miranda, Opher 
& Opher (1998) concluded that field seeds as high as B 4.5 x 10"^° G can be created with 
coherence scales of order of 100 kpc. However, their mechanism requires efficient SN shell 
fragmentation, which is not expected to occur in such low metallicity environments (Fcrrara 1998). 
Furthermore, most of their magnetic flux is produced by objects with mass W^Mq forming at 
redshift z 300, which is clearly at odds with all currently viable cosmological models. 

We should also mention that battery mechanisms could have operated on much smaller scales, 
such as in accretion disks around massive compact objects, and seeded young galaxies and the 
IGM through outflows and jets (Daly & Loeb 1990). 

Observations provide little guidance in choosing among these theoretical models. Attempts 
to determine intergalactic magnetic fields have yielded uncertain results, or at best upper limits. 
A very clean determination of primordial i3-ficlds could be accomplished through measurement 
of Faraday rotation in the polarization of the CMB (Kosowsky &; Loeb 1996). This will become 
possible with forthcoming experiments. 

An experimental upper limit Bq ^ 10^^ {Qigm/ 0-01) G (see review by Vallcc 1997 and 
references therein) is found on the intergalactic B-&cld averaged over a length scale comparable to 
the horizon. This upper bound would increase by a factor of 3 if the field were coherent on 10 Mpc 
scales (Kronberg 1994). 

Finally, there is some evidence for ^uG magnetic fields in young galaxies, as traced by 
high redshift (up to z ^ 3) quasars (Kronberg, Perry Sz Zukowsky 1990) and damped Lya 
systems (Wolfe, Lanzetta & Oren 1992). This value might reflect the saturated level dictated by 
equipartition with turbulent energy (Field 1995), but it requires a rather high initial seed field 
value in order for the dynamo to amplify the field up to that strength. At that epoch, a dynamo 
working in a galaxy like the Milky Way should be provided a seed B ^ 10~^^ — 10~^^ G (Beck et 
al. 1996). 

The purpose of this paper is to explore magnetogenesis in cosmological ionization fronts 
using new simulations of the reionization of the universe by stars in protogalaxies. Our results 
include, in a self consistent manner, the wealth of density structures present at the epoch of 
galaxy formation (at least down to a scale imposed by numerical resolution). We corroborate the 
results of SNC to order of magnitude, place their suggestion on a firmer ground, isolate the most 
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Table 1. Numerical Parameters 



Run 


N 


Box size 


Baryonic mass res. 


Spatial res. 


A 


128^ 


4/i"^ Mpc 


IO^-'^Mq 


l.Oh-^ kpc 


B 






1O^-'^M0 


1.5/i"ikpc 


C 


128^ 


2/1-1 Mpc 


10^-^ Mq 


0.5/1-1 j^pp 



important mechanisms, and calculate the structure of the magnetic field in greater detail. 

In §2 we describe the simulations, §3 gives the results, and §4 is a summary of the calculations 
and their limitations, as well as a discussion of their astrophysical relevance. 

2. Method 
2.1. Simulations 

We use cosmological hydrodynamic simulations performed with the "Softened Lagrangian 
Hydrodynamics" (SLH) code and reported in Gnedin (2000). These simulations include 3D 
radiative transfer (in an approximate implementation) and are therefore suitable for following 
magnetic field generation. In addition, the simulations include a phenomenological description of 
the star formation process, based on the Schmidt law and on a full time-dependent treatment of 
atomic and molecular physics in a plasma with primeval composition. 

We adopt a currently fashionable variant of the CDM+A cosmological model with f^o = 0.3, 
0,\ = 0.7, h = 0.7, and ^Ib = 0.04. We have run three simulations with varying resolution, to 
estimate the level of numerical convergence. The numerical parameters of all three runs are listed 
in Table ||. Runs A and B are stopped at z = 4 because at lower redshift the simulation box cannot 
be considered as a representative region of the universe: in fact, the rms density perturbation 
on the scale of the box size exceeds 0.25. Run C was stopped at z = 6.5 because it is used for 
estimating numerical convergence only. 

We will use run A as our reference run, whereas runs B and C will serve as a measure of the 
theoretical uncertainty of our calculations. 

The placement of the simulation box of run A is designed in such a way as to marginally 
resolve the characteristic filtering scale below which the baryonic perturbations are smoothed due 
to the finite gas pressure. Thus, run A includes essentially all (about 93%) of the small scale power 
that is present in the initial conditions (Gnedin 2000). Since the small scale perturbations can 
contribute significantly to the total production rate of the magnetic field, we thus ensure that we do 
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not miss a substantial portion of the total production in the regions with moderate overdensities. 
However, inside the virialized objects small scale power can be generated gravitationally, and since 
this power is not resolved, we are almost certainly underestimating the true value of the magnetic 
field. We argue below, however, that the magnetic field on large scales is generated primarily by 
structure on large scales. Thus, our calculation should produce a good approximation to the large 
scale fields produced by magnetogenesis in ionization fronts. The next stage, amplification of these 
primordial fields by dynamo action inside virialized objects, is outside the scope of this study. 



2.2. Magnetic Field Equation 

The evolution equation for the magnetic field with the Biermann battery source term can be 
written as follows: 

^ + {v-V)B = -B divv + {B ■V)v + —VT X Vn^ + —Ue^V x v (1) 
at rieC e 

where rig is the electron number density, T is gas temperature, and 

with U-y being the energy density of the CMB, and cjt the Thompson cross section. 

The left hand side of equation (||) describes a simple advection of the magnetic field together 
with the fluid. The two first terms on the right hand side of the field (induction) equation can be 
rewritten in index notation in the following form: 

-B^ + B^v]^ = -\&v\^ + B\v\^ - i^.^) = + B^a)^ 

where Dp/dt is the Lagrangian derivative of the gas density, and is the traceless part of the 
velocity tensor. These two terms describe two physically different processes of amplification of the 
magnetic field : compression and stretching respectively. We discuss the importance of those two 
terms below. The third term in equation is the Biermann battery, and the last term is due to 
Compton drag of CMB photons on the free electrons. We ignore the Hall effect and consider ions 
and neutrals moving with the same fluid velocity v. 

From general physical considerations, we expect that the battery term always dominates in 
ionization fronts, and is small everywhere else. It is useful to estimate the magnetic field produced 
by the battery term as an ionization front sweeps through the gas at speed Vf. Integrating 
equation (||) in time, keeping only the battery term and assuming that Vng is in the direction of 
propagation of the front while VT is perpendicular to Vng and is produced by some large scale 
inhomogeneity, we find 

B ~ — , 

evflg 
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where Ig is gradient length scale for T. In order to get some physical feeling for this field strength, 
we introduce the ionization time ti = Ig/vf, the ion gyrofrequency = eB/niiC, and the ion 
thermal velocity vth = {k^T /niiY/'^ . In terms of these quantities, 



Or, in terms of the thermal ion gyroradius rj = Vth/'^d, 

Ti Ig ~ — . 

Vth 

In contrast to the usual astrophysical situation, in which uj~^ and are much smaller than any 
global time scales or length scales, the battery produces a field so weak that the magnetization 
scales and global scales are comparable. But note that small scale structure (small Ig) produces 
stronger fields. 

The ratio of the photon drag term to the induction term, assuming that the vorticity is 
appreciable, is 

~ ^, (2) 

where uJce ~ 10'''-Bs~^ is the electron gyrofrequency. Since = 5 x 10~^°(1 + z)'^s~^, the 
induction dominates photon drag for 

S > 5 X 10-27(1 + ^)4 0. (3) 



Using the outputs from the simulations, we solved equation (||) numerically using clouds- 
in-cells (Efstathiou et al. 1985) to advect the right hand side of equation (||) over the SLH 
quasi-Lagrangian mesh, and implementing the Lax averaging for the stretching term to ensure the 
stability of the scheme. We chose the time-step so as to ensure numerical convergence with respect 
to time discretization. 

We ignore the back reaction of the magnetic field because the magnetic pressure for the field 
strength we obtain [B ~ 10^^^ G, B'^/8it ~ 10^'*'^ erg/ cm'^) is far smaller than the gas pressure 
even in the low density IGM {P ~ 10^"'^^ erg/ cm^ at z ~ 7). 



3. Results 

3.1. General Outlook 

We can now look at general features of the magnetic field generated during the reionization. 
Figure |l| shows the evolution of the mass and volume weighted mean magnetic field strength as 
well as the mean free path to ionizing radiation for all three simulations discussed in this paper. 
The difference between the curves illustrates the level of numerical convergence achieved in the 
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Fig. 1. — Evolution of the mass {bold lines) and volume {thin lines) weighted mean magnetic field strength 
{top panel) and the comoving mean free path to ionizing radiation {bottom panel) for runs A {solid lines), 
B {dotted lines), and C {dashed lines). 

simulations, and can be considered the theoretical error of our simulations. In particular, the 
differences at high redshift are due to the fact that our reference run A simulation still misses 
some 7% of the small scale power, which translates into a ~ 15% shift in redshift. However, the 
saturation levels of all three runs are about the same. Since run C has a twice smaller cell size 
(and, therefore, approximately four times less numerical diffusion), the agreement in the saturated 
level of the produced magnetic field demonstrates that numerical diffusion does not lead to a 
substantial underestimate of the magnetic field strength produced in our simulations. This test 
is not however completely conclusive because of the excess small scale power present in run C, 
and therefore we cannot exclude the possibility that numerical diffusion affects our results on a 
several tens of percent level. However, since we focus in this paper on the general semi-qualitative 
description of the generation of the primordial magnetic field, we can tolerate a factor of two 
uncertainty in our calculations. 
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Fig. 2. — Evolution of the mass (hold lines) and volume (thin lines) weighted mean magnetic field strength in 
the full calculation (solid lines), with the Compton drag term omitted (long-dashed line), with the stretching 
and the Compton drag terms omitted (short-dashed line), and with the compression and the Compton drag 
terms omitted (dotted line). 

We point out here that the mean field strength generated in our simulations exceeds the 
value found by Kulsrud et al. (1997) by about two orders of magnitude. We will elaborate on this 
difference in the next subsection. 

The evolution of the mean field clearly exhibits two distinct regimes: before the overlap of the 
H II regions at z = 7 (characterized by a sharp increase in the mean free path of ionizing photons) 
the mean field grows with time, but after the overlap its growth slows down appreciably. We 
attempt to understand this behavior by considering which terms in equation (0) are dominant at 
a given time. Figure || shows four solutions to equation ([l|) for the large simulation: the solid lines 
mark the full solutions (either mass- or volume- weighted) , the same one as shown in the previous 
Figure; the long-dashed line shows a solution with the Compton drag term omitted; the dashed 
line shows a solution with both the Compton drag term and the stretching term omitted; and the 
dotted line shows a solution with both the Compton drag term and the compression term omitted. 
We can immediately conclude that stretching by itself makes an insignificant contribution, whereas 
compression is the dominant term after the overlap of the Hii regions, and the Biermann battery 
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Fig. 3. — The joint mass-weighted distribution of the gas density and the magnetic field strength at z = 4. 

is responsible for the initial growth of the field during the pre-reionization stage. The Compton 
drag term is only important for z > 16, when the magnetic field is less than about 5 x 10"^'^ G. 
Using the estimate we find that at this moment the induction term is about 10% of the 
Compton drag, and the comparison between the dotted and solid lines in Fig. ^ shows that the 
Biermann battery at that time is on average about 10 times more important than the induction 
term, in agreement with our estimate @. 

Since compression is the dominant term at lower redshifts, we expect that the magnetic field 
has to be closely related to the gas density. The battery term has a similar effect, as ionization 
fronts tend to move slowly in regions of high gas density. Figure ^ illustrates this point. In 
this Figure we show the joint mass-weighted distribution of the gas density and the magnetic 
field strength for all fluid elements in the large simulation. As one can see, there is a strong 
correlation between the gas density and the magnetic field, albeit with a considerable (two orders 
of magnitude in magnetic field strength) scatter. The most probable value for the magnetic field 
strength (the peak of the distribution) rises with density slightly more steeply than the first power. 
We will give an explanation for this feature in the next subsection, where we discuss the main 
mechanism for generating the primordial magnetic field. 

In Figure Q we show a thin slice through the simulation box taken at a random place. In four 
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Fig. 4. — A thin slice through the simulation volume at z = 4, showing in four panels the logarithm of 
the neutral hydrogen (upper-left) and the gas density (lower-left) (in units of the average density of the 
universe), the gas temperature (lower-right) , and the magnetic field strength (upper-right) . 

panels we show the gas density, temperature, neutral fraction, and the magnetic field strength 
(all stretched logarithmically) .0 Again, visual inspection demonstrates that the magnetic field is 
strongly correlated with the gas density. It is also correlated with the gas neutral fraction and 
the temperature, as these quantities also correlate with the gas density in the moderate density 
regime (cosmic overdensity 6 ^ 10). 

Finally, in Figure |5| we show the magnetic field strength again (just as in Fig. ^) together with 
the three cartesian components of the magnetic field. Since the components can be both positive 
and negative, we use the following stretch to map the value of the component in the range from 
— 10~^^ G to +10~^^ G to a quantity in the range from —3 to +3: 

stretch = sign(Si) log(l + 10^^ G"^[Si|). 



As one can see, the magnetic field is highly ordered on scales of the order of 1 Mpc, and magnetic 
structure is highly correlated with density structure. As pointed out above, this is the result of 



^Additional visualization of the results of the simulations including MPEG movies can be found at the following 
URL: http : / / casa. Colorado . edu/'gnedin/GALLERY/magf i_p .html. 
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Fig. 5. — A thin slice through the simulation volume at z = 4, showing the magnetic field strength in the 
upper-right panel (just as in Fig. ^) and three components of the magnetic field in the other three panel. 

the compression occurring as cosmic structures form. As we have already discussed, this is not 
an artifact of finite numerical resolution, since our simulation resolves all the relevant scales in 
the moderate density regime. Only inside the virialized clumps do we encounter scales which we 
are unable to resolve, and there our calculation of the magnetic field strength becomes a lower 
limit, on top of which small scale dynamo processes will generate additional field on galactic 
and subgalactic scales. Since these latter processes cannot be simulated reliably with the current 
resolution of our simulations, we cannot consider them in this paper. On the other hand, the 
extragalactic magnetic fields are calculated reliably in our simulations for a given cosmological 
model. 



3.2. Main Mechanisms for Generating the Magnetic Field 

Since the Biermann battery is most efficient when the gradient of the electron density is 
perpendicular to the temperature gradient, we can identify two main mechanisms where the 
battery is most efficient. But before we can do so, let us briefly describe the main stages of 
reionization (Gnedin 2000). We also point out that the reionization process is affected by several 
feedback mechanisms which are discussed in Ciardi et al. (2000). 
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The reionization starts with ionization fronts propagating from proto-galaxies located in high 
density regions into the voids, leaving the high density outskirts of the protogalaxies still neutral, 
because at high density the recombination time is very short, and there are not enough photons to 
ionize the high density regions. As an ionization front expands into previously neutral material, 
it leaves behind high density regions whose ionization requires more photons than currently 
available. This stage of the reionization process can be called "pre-overlap" , and it extends over a 
considerable range of redshifts Az ~ 5 around z ~ 10. During this time the high density regions 
close to the source slowly become ionized, whereas more distant high density regions remain 
neutral. 

By z ~ 7 (for the cosmological model under consideration) the Hii regions start to overlap. 
As a result, an average place in the universe can see an increasing number of sources, and the 
ionizing intensity starts to rise rapidly. The process of reionization enters its second stage, the 
"overlap", which is quite rapid (Az ~ 1). As the ionizing intensity is rapidly increasing, the last 
remains of the neutral low density IGM are quickly eliminated, the mean free path increases by 
some two orders of magnitude over a Hubble time or so (see Fig. 1), and voids become highly 
ionized (neutral fraction of the order of 10~^). The high density regions at this moment are still 
neutral, as the number of ionizing photons available is not sufficient to ionize them. 

After the overlap is complete, the universe is left with both highly ionized low density regions 
and some of the high density regions (which happened to lie close to the source, where the local 
value of the ionizing intensity is higher than the background). High density regions far from any 
source remain neutral. This stage can be called "post-overlap". As time goes on, and more and 
more ionizing photons are emitted, the high density regions are gradually eaten away. 

Figure |^ shows now a cartoon version of the first out of two most important mechanisms for 
generating the primordial magnetic field. The mechanism takes place when the first I- front breaks 
through the protogalaxy. This happens first in a narrow range of angles, where the gas density 
happens to be the lowest. As the I-front propagates into the low density IGM, the gradient of the 
electron density (the surface of the I-front) becomes approximately orthogonal to the temperature 
gradient, and the battery becomes maximally efficient. 

Figure |^ illustrates this mechanism with an example from the simulation. Each of the frames 
is similar to Fig. ^, except that it is only one fifth the box size (0.8/i~^ Mpc) on a side. Two 
panels show a protogalaxy just before and after the I-front broke through the gaseous halo.0 
The magnetic field has increased by almost an order of magnitude during the short time interval 
(Af = 75 Myr) between the two snapshots. 

The magnetic field produced by this mechanism has toroidal topology, wrapping around 



^The slice through the simulation box is chosen randomly, and is not passing through the center of the object. We 
are therefore seeing a cross section of Fig. ^ taken off center and at an angle with respect to the axis of the ionized 
cone. 



- 13 - 



Electron 


density gradient 




MPI 1 t Vri 1 






loni'/ed 










Temperature 


gradient 


Stars 



Fig. 6. — A cartoon, illustrating the mechanism for generating the primordial magnetic field during the 
break through of the I-front from the proto-galaxy before the overlap of the Hii regions. Enclosed contours 
show the regions of progressively higher density. 

the H II region. This can be seen by inspection of the battery term: Vng is directed across the 
ionization front, while VT follows the large scale stratification of the protogalaxy and is directed 
along the major axis of the cone shaped H II region. The magnetic field is oriented along the cross 
product of these gradients, i.e. the field encircles the H II region. The reversal of the out-of-plane 
component of the field across the H II region is indeed observed in the simulation. 

The second mechanism is important when the I-front propagates across the high density 
neutral filament. This is shown in Figure ^ as a cartoon version and further illustrated in Figure ^ 
with two snapshots from the simulation. Both, the cartoon and the snapshots, demonstrate how 
this mechanism operates in the post-overlap stage, when the low density IGM is already ionized. 
However, it is clear from the cartoon that the ionization state of the surrounding low density IGM 
is irrelevant, and that this mechanism works in the pre-overlap stage as well, when the filament 
is embedded in the neutral rather than ionized gas. In general, however, I-fronts propagate much 
more slowly across the high density filaments than in the low density IGM, so it is likely that the 
IGM around the filament is already ionized when the I-front propagates across it. 

It is important to keep in mind that the filaments are neutral not because they are 
self-shielded, but because they are dense enough for the recombination time to be shorter than the 
photoionization time, i.e. even if filaments are neutral, they are nevertheless photoheated. Since in 
the IGM with moderate overdensities {6 ^ 10) there exists a tight relation between the gas density 
and the gas temperature (the so called "effective equation of state", Hui & Gnedin 1998), even a 
filament with quite moderate overdensity has a temperature gradient pointed outward, the most 
efficient orientation for Biermann battery. 
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Fig. 7. — A thin slice through the simulation volume similar to Fig. ^ but showing l/5th of the box 
(O.Sft.^^ Mpc on a side), aimed at illustrating the generation of the magnetic field during the break-through 
of the I-front. Two redshifts are shown: z = 11.5 (a) and z = 10.1 (6). The region around the I-front is 
emphasized by the white circle. 



Ionized 




Temperature 
gradient 



Electron 
density 
^gradient. 



Jonization 
front 



Fig. 8. — A cartoon, illustrating the second mechanism for generating the primordial magnetic field during 
reionization: the I-front crossing a neutral high density filament after the epoch of overlap. Enclosed contours 
show the regions of progressively higher density. 
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Fig. 9. — A thin slice through the simulation volume similar to Fig. ^ but showing l/25th of the box 
{0.8h^^ Mpc on a side), illustrating the generation of the magnetic field during the propagation of the I-front 
across the high density filament. Tho redshifts are shown: z = 5.9 (a) and z = 5.5 (&). The region around 
the I-front is emphasized by the white circle. 

The first mechanism operates during the formation of a virialized objects, i.e. in higher 
density environments than the second mechanism. This fact explains why the field increases with 
density faster than p^^^, which would be the predicted law if the simple compression were the 
dominant term in equation (||). Since high density regions collapse first, the field per unit mass 
inside the protogalaxies is larger than the field per unit mass in the lower density filaments and 
the IGM. 

The existence of these two mechanisms is also the source of our disagreement with Kulsrud et 
al. (1997), who found fields which are two orders of magnitude less than the value we find. Since 
Kulsrud et al. (1997) did not include radiative transfer in their simulations, they were not able 
to follow the I-fronts, and therefore completely missed these two important mechanisms for field 
generation. Their field was therefore generated by chance misalignments between the electron 
density and temperature gradients. These misalignments are also present in our simulations, but 
they can be a dominant source of the magnetic field only in the low density IGM, where none of 
the mechanisms described above operates. Thus, in the low density IGM we should be finding 
values for the magnetic field strength similar to the ones found by Kulsrud et al. (1997), i.e of the 
order of 10~^^ G, which is indeed the case, as can be seen from Fig. |3|. 
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3.3. Primordial Field in the Protogalaxies 

As we have mentioned above, our simulations do not have enough resolution to model the 
interstellar medium, and therefore cannot describe the field generation in the galactic ISM. Thus, 
the values of the magnetic field strength that we find inside the protogalaxies should only be 
considered as the truly primordial field, generated during the period of reionization. The true field 
inside the protogalaxies at, say, z = 4, could be significantly higher than the values we obtain (of 
the order of 10~^^ G), although it would be coherent only on small scales. Here, we have in mind 
small scale fields generated directly by the battery mechanism, as ionization fronts propagate 
through the clumpy molecular clouds in which stars are born. However, these small scale fields do 
not contribute as much to the magnetization of protogalaxies as large scale fields do, as we show 
by the following rough argument. 

Suppose that the battery acts on small clouds or filaments of characteristic size Do and 
density ric, and produces a field B{Do). If most of the mass in the star-forming complex 
is in these small clumps, then when the clumps eventually disperse their mean size will be 
D ~ Do{nc/ (n))^/^ , where (n) is the mean interstellar density. Because of conservation of 
magnetic flux, the magnetic field in the dispersed clumps is B(D) ~ B{Do){{n) /nc)"^^^ . Space 
is then divided into independent magnetic domains of size D, and the fields within the domains 
are randomly oriented with respect to one another. Over a region of size I > D, sl random walk 
argument shows that B{1) B{D){D /Vf/'^ (Hogan 1983). In terms of the original battery field, 
B{1) ^ B{Do){{n)/n,fl\Do/lfl\ 

We now compare B{1) computed as the sum of fields in small domains with the coherent 
field Bc{l) which is generated by propagation of a single I- front. From the scaling of the battery 
with system size and I-front velocity vj we expect Bc{l) ~ B{Do){Do/l)[v f{Do) /v f{l)]. The ratio 
BS)/B{1) is then Bc{l)/B{l) ~ {1/ Dof/'^{nJ {n)f/^[vf{Do) /vf{l)]. This expression shows that 
unless the ionization fronts in the clumps propagate much more slowly than the fronts on large 
scales, the coherent field dominates the superimposed incoherent fields. Although the clumps can 
be very dense, which slows down the fronts, they are also much closer to the source of ionization, 
which speeds them up. Therefore, the ionization front battery on large scales is more important 
for the generation of a galactic field than the battery on the scale of molecular clumps in star 
forming regions. However, if the locally generated battery field is amplified sufficiently rapidly to 
equipartition with the turbulence, it could play a role in subsequent generations of star formation 
in the protogalaxy. 

Since we follow the intergalactic field with sufficient precision, we can separate the whole 
problem of the generation of the magnetic field into two distinct phases: the formation of the 
primordial field during reionization, which we can reproduce with our simulations, and the dynamo 

action inside the protogalaxies, which wc cannot resolve. The latter process, though, is restricted 
to the protogalaxies only, and therefore could be, at least to a first approximation, considered 
independently from the rest of the universe. 
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Fig. 10. — Evolution of the stellar mass (upper-left), the gas mass (upper-right), the rms magnetic field 
strength (lower-left), and the absolute value of the mean z-component of the magnetic field of the most 
massive virialized object within the 4/i~^ Mpc simulation box (solid lines), the 10th most massive object 
(long-dashed lines), the 50th most massive object (dot-dashed lines), and the 250th most massive object 
(dotted lines). 

Figure |l^ shows the time evolution of the stellar mass, the gas mass, the rms magnetic 
field strength and the absolute value of the mean z-component of the magnetic field (x- and y- 
components look very similar) for four representative objects from our simulation. We notice that 
the rms field traces the evolution of the mean cosmic field from Fig. |l|, and saturates at the values 
of about (1 - 2) X 10"^^ G in the most massive objects. The absolute value of the z-component is 
about 10 times smaller, implying that only 10% of the total field is in an ordered state over the 
scale of an object^. 



Finally, in Figure 11 we show the final (i.e. saturated) value of the rms magnetic field 



strength versus mass for all virialized objects from our simulation. The most massive objects with 



®The dips in the absolute value of the ^-component are due to rotation of an object, when passes through zero 
and changes sign. 
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Fig. 11. — The final rms magnetic field strength versus baryonic mass for all virialized objects from our 
simulation. 

baryonic mass in the range 10^ — 10^° Mq all have rms magnetic field strength about 10~^^ G with 
a factor of a few scatter. One can also notice a correlation between the baryonic mass and the 
rms magnetic field strength, albeit the scatter is quite large. This correlation occurs for the same 
reason as the steep B vs. p dependence: the 1-fronts are more efficient in producing magnetic 
fields while they are breaking through the neutral halos of the protogalaxies than when they are 
propagating across filaments. 



4. Conclusions 

We have studied the generation of magnetic field during the reionization of the universe 
outlining the main mechanisms, arising in ionization fronts, responsible for magnetogenesis. Our 

results allow us to conclude that the cosmic seed magnetic field closely traces the gas density, and 
that it is highly ordered on megaparsec scales, with a mean strength Bq 10~^^ G. We find that 
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in protogalaxies the field is an order of magnitude or more larger than the mean field; this is an 
upper limit because unresolved small scale structure in the protogalaxies could generate stronger 
fields, and amplify them rapidly. 

The magnetogenesis mechanism proposed here owes its robustness to both the relatively 
simple physics involved {i.e. ionization fronts) and the strong experimental evidence that 
reionization of the universe has indeed occurred. Although we have assumed here that reionization 
is driven by stellar sources, which appears to be consistent with available experimental data on the 
metallicity of the IGM (Gnedin &; Ostriker 1997; Ferrara, Pettini &; Shchekinov 2000), quasars - if 
present prior to reionization - are likely to produce a similar effect as they ionize the surrounding 
intergalactic medium. 

This field strength is about 10 orders of magnitude too small to have a dynamical effect on 

galaxy formation. One way to see this is to recognize that the Alfven speed in a medium with 
n = 10^^ cm^^, B = 10^^^ G is 2 x 10^^ cm/ s. This field is also too small to have much effect on 
microscopic transport processes such as thermal conductivity. The relevant parameter here is uJcTc, 
where lOc and Tc are the gyrofrequency and Coulomb collision time, respectively, and their product 
must exceed unity for thermal conduction to be suppressed. In a plasma with n = 10~^cm~^, 
T = 10^ K, B = 10~^^ G, this parameter is about 3 x 10~^ for electrons. Similarly, the gyroradius 
of a high energy cosmic ray propagating through this field is extremely large: a 10^^ eV proton in 
a 10^^^ G magnetic field has a gyroradius of 10^^ Mpc, so high energy intergalactic cosmic rays 
would be essentially unmagnetized. 

However, such a field, even reduced by (1 + z)^ from its value at formation, might be 
detectable through its effects on the arrival times of 7-rays from extragalactic sources (Plaga 1995, 
Lee, Olinto, & Sigl 1995, Waxman &; Coppi 1996). The effect of the field on the polarization of 
the CMB would be much less than that modeled by Kosowsky Sz Loeb (1996), who required fields 
9-10 orders of magnitude larger to find a detectable magnetic signature. 

This field could explain the strength currently observed in the Milky Way and in external 
galaxies, provided it is further amplified by a some other process, as for example a dynamo which 
has exponentiated 30 times. The protogalactic seed fields are just barely large enough to have 
been amplified by a dynamo to /iG strengths by z ~ 1, as suggested by the observations of damped 
Lya systems. 
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